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Abstract 



We discuss some puzzling aspects of the narrow hidden charm resonance 
that was recently discovered by the Belle Collaboration at mass 3872 MeV. 
In order to determine its quantum numbers, a crucial piece of information is 
the spin of the dipion in the decay final state -k+it^ J/ip. We give the angular 
distributions and correlations of the final particles in the decay which will 
provide this information about the nature of this resonance. 
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I. INTRODUCTION 



A new narrow resonance has been discovered at mass 3872 MeV by the Belle Collabora- 
tion [1] through the process, 

B ± -> K ± X(3872) -> K^+ir-J/ip. (1) 

Following Belle [2], we denote this resonance tentatively by X(3872) in this paper. The 
invariant mass of 7r + 7r~ extends to the upper end of its kinematical boundary (~ 775 MeV), 
but it is not known at present whether the dipion is in s-wave, p-wave, or even ci-wave. 
Experiment does not exclude the possibility that the dipion in the p mass region is actually 
in s-wave. We should keep in mind that the s-wave nn scattering cross section rises rapidly 
starting just below the p mass. In fact, this partial-wave or spin of the dipion provides the 
most important clue to the quantum numbers of this resonance. Search of the radiative 
decays X(3872) — > jxa has so far not produced a positive result [2]. As the upper bound 
set on the raditative decay branching becomes more severe, it will impose a strong constraint 
on the nature of the resonance. 

The most likely candidate for X(3872) is an excited charmoniun state. In fact, the ex- 
perimental study of the 7i + 7i~ J/ip final state in B decay was motivated with search of the 
excited charmonia [3]. Meanwhile, closeness of the DD* threshold to 3872 MeV suggests an- 
other explanation, a loosely bound molecular state of DD and DD* [4]. However, both the 
charmonium and the molecule interpretation may encounter some problems. In this paper 
we are not rigidly constrained by the potential model predictions on charmonia since un- 
certainties are large for the excited charmonium states near the open charm thresholds. We 
will be open minded about dynamics of the molecules. Instead we narrow down the possible 
quantum numbers of X(3872) with the information extracted from the Belle experiment and 
then focus on possibility of determining the quantum numbers purely experimentally by the 
decay angular distributions and correlations. 



II. CURRENT EXPERIMENTAL INFORMATION 

The Belle experiment [1] has so far imposed the following constraints on this resonance. 
The width is narrow (< 2.3 MeV [2]) despite the ample phase space (p cm — 500 MeV) 
for the decay X(3872) -> DD. For comparison, the width of ^"(3770) -> DD is 24 MeV 
with p cm = 242 MeV. Provided that the decay into DD is forbidden by selection rules of 
quantum numbers rather than by some unknown dynamical suppression, we expect that 
X(3872) should have an unnatural spin-parity, 

J P = 0-,l + ,2----, (2) 

or an unnatural spin-charge-parity, 

j pc = 0+ - r + 2 +-..._ (3) 

To select a right J PC out of these choices, the dipion quantum numbers provide the most 
important clue. If n + n~ forms a scalar dipion of J^p = ++ or a tensor dipion of J^p = 2 ++ , 
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charge parity C of 7r + 7r~ J/?/> is negative and isospin is / = or 2. We do not consider 1 = 2 
further in view of lack of any candidate. Combining 7r + 7r~ with J/ip in relative orbital 
angular momentum L, one can make the unnatural spin parity and charge-parity states 
of + ~, 2 ±_ • ■ • with a scalar or a tensor dipion. If ir + ir~ forms a vector dipion of 

= 1 , charge parity of tt + tt^J/^ is positive and isospin is / = I. The most relevant 
unnatural spin-parity state is 1 ++ with L = in this case. These unnatural quantum number 
states are listed in Table I. Since only a limited phase space is available for the J /ip 
decay, the Table includes only the cases of L < 2 for the scalar dipion, and L < 1 for the 
vector and tensor dipion. 



III. CHARMONIA 

Since / = for the charmonia, the folllowing quantum numbers are selected for charmo- 
nium candidates: 

t pc_ f2-- (L = 0,2) m 
" \ (L = 1). W 

The charmonia carrying these quantum numbers are 3 D 2 (2 ) and 1 Pi(l H ). 

The mass spectrum calculations by potential models [5,6] suggest that the 3 D 2 state 
should be much closer to the l 3 Di state, which is believed to be -?/>"(3770). The mass 
of the 1 3 D2 state was predicted at 30~60 MeV lower than 3872 MeV [5]. However, the 
charmonium levels near the DD and DD thresholds are subject to large uncertainties due 
to the open charm channel coupling and the large relativistic corrections (y^ {v 2 /c 2 ) ~ 0.54 
in one estimate [6]). Therefore we are not too concerned with the discrepancy between 
experiment and the potential model expectation for the masses in the case of those higher 
excited charmonia. 

The radiative decay into 7Xcj is allowed for both 2 and l +_ by quantum numbers. The 
decay 3 D2(2 ) — > 7Xcj is the El transitions and should occur without suppression. Let us 
compare it with a similar radiative decay of ip{2S). We know Y{ip{2S) — > xa) — 26 keV and 
T(i/j(2S) — > n+iT^J/'ip) ~ 93 MeV from experiment. The strong decay i>{2S) — > 7r + 7r - J/ip 
occurs into a scalar dipion with s-wave (L = 0) relative to J/ip. Apart from difference 
in the dipole strength (r) between ip(2S) — > Xcj an d l 3 D 2 — > Xcj, the spin factors and 
phase space corrections combined enhances r( 3 D2(3872) — > 7Xci) by a factor of ~ 5 relative 
to r(ip(2S) — > 7Xci)- The dipole matrix element is sensitive to the coupling to the p- 
wave DD channel. Comparison of the strong decay is at least as sensitive or even more 
uncertain since the decay X(3872) — > 7i + 7i^J/ip involves d wave either in the dipion or in the 
relative orbital angular momentum L. The rescaling factor of the three-body phase space 
can be as large as ~ 10, which is severely compensated by the rf-wave suppresion factor, 
| Pttt,- / -Etttt 1 4 f° r a scalar dipion and Ip^/i^ 4 for a tensor dipion. It depends sensitively on 
the dipion mass distribution. Despite the larger phase space for X(3872) — > tt + tt^J/^ 
than for tp(2S) -> tt+tt-J/V', the net result is most likely that T(X(3872) -> tt+tt-J/^) 
is smaller than r(ip(2S) — > 7r + 7r~ J/ip). With this order of magnitude consideration we 
expect J B(X(3872) -> -f X ci)/B(X(3872) -> tt+tt-J/^) > 1. Then the non-observation of 
the radiative decay X(3872) — > ^Xci [2], 
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5(X(3872) -> 7x c i)/5(X(3872) -> Tr+Tr" J /if;) < 0.89 (90% C.L.), ( 5 ) 

poses a problem on the 3 D 2 assignment. The Belle Collaboration [2] compares the upper 
bound of Eq. (5) above with the potential model prediction of 5 [3] . Considering the large 
uncertainties involved in the potential model calculations, however, it is prudent to take a 
"wait and see" attitude. 

The other possibility of the radially excited 1 Pi(l + ~) encounters a larger discrepancy 
with the mass prediction of the potential models. The mass 3872 MeV is roughly 100 MeV 
lower than the potential model. Nonetheless we do not reject 2 1 P 1 at present for the same 
reason stated above for 1 3 D 2 . The radiative decay 2 l P x — > 7Xcj is an Ml transition that is 
caused by spin flip. Since the spatial wave functions are orthogonal between x c j(1 3 Pj) and 
2 l Pi in the nonerelativistic limit, the Ml transition is considerably weaker than the allowed 
El transitions. The strong decay 2 1 P± — > it + it~ J/ip occurs into a scalar dipion with L — 1 
instead of L = 2 for 1 3 D 2 — > it + it~ J/ip. Therefore the radiative branching fraction is much 
smaller, relative to the it + it~ J/ip decay branching, for 2 l P\ than for 1 3 £>2- It can be easily 
consistent with the current upper bound. If the non-observation of the radiative decays into 
IXcj continues to a higher precision, we should consider 2 X P\ as a better candidate than 
l 3 .D2- We should keep in mind, however, that even the lowest l^Pi charmonium h c has not 
yet been reported in B decay despite continuing searches by both the Belle and the BaBar 
Collaboration. 

To wit, if X(3872) is a charmonium, it should be either the 1 3 -D 2 (2 ) state or the 
2 1 P 1 (1+-) state. 



IV. MOLECULES 

The idea of loosely bound molecule states of two hadrons has been entertained for a 
long time [4], but no meson has so far been positively identified as such a state. The close 
proximity of the mass 3872 MeV to the DD* thresholds (3971.2 ± 0.7 MeV for D°D*° and 
3979.3 ± 0.7 MeV for D + D*~) has prompted reconsideration of this possibility for X(3872) 
[7]: 

(1/V2)(DD* ±DD*). (6) 

Hereafter we refer to the charge-parity eigenstates of the molecule states simply by DD* . 
Then the following unnatural spin-parity and charge-parity states can be formed as a 
molecule: 

T pc _ J l +± (W = 0), 



-\o-M-+,2-± = [7) 

Since charge parity of 7r + n~J/ilj is correlated to isospin by C = — (— l) 7 , the molecule should 
be in / = for C = — and in / = 1 for C — +. For the relative orbital angular moementum 
of DD , the closeness of 3872 MeV to the DD threshold strongly favors L D p* = over 
L D jf = 1. A qualitative argument of dynamics based on the one-pion exchange favors 
formation of the 1 = molecules over the 1 = 1 ones, in particular, _+ and 1 ++ [8]. 
However, the states of 0~ + and 1 ++ with 1 = are not found in Table I since they are 
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inconsistent with the decay into 7r + 7r~ J '/ip. Is it possible that the states of 0~ + and 1 ++ 
with 1 = decay into the 1 = 1 channels of 7r + 7r~ J/ip with isospin symmetry breaking. 
The answer is affirmative, but such a state cannot feed a scalar dipion in n + n^J/ip because 
of C invariance. In this case, the dipions observed with mass below the p region would be 
entirely due to the vector dipion produced by isospin breaking. 

We do expect a large isospin violation in a molecular bound state of DD* at 3872 MeV, 
since the mass 3872 MeV almost coincides with the D°D threshold but as much as 7 
MeV below the D + D threshold. No matter how large isospin violation is, however, no 
single resonance can produce both a scalar dipion and a vector dipion in TT + n" J/ip by C 
invariance: Any strong interaction resonance of zero net flavor must be an eigenstate of 
charge conjugation even when isospin is broken. Since the states of / = and 1=1 have 
opposite C, a single resonance cannot be responsible for both the 1 = final state (a scalar 
dipion) and the 1 = 1 final state (a vector dipion) of 7r + 7r~ J/ip. In order to feed both 
the scalar/tensor and the vector dipion, there must be two (almost perfectly degenerate) 
resonances of opposite charge parities, one with 1 = and odd C and the other with 1=1 
and even C. Forming two such molecules in degeneracy in hidden flavor channels would be 
a highly unlikely accident. With these observations, we should consider only a molecule of 
either J PC = l + ~ with I = or J PC = 1 ++ with I — 1 as the molecule candidates among 
the entries in Eq. (7). 

It is likely that the branching fractions of the radiative transitions to 7Xcj are naturally 
small for the DD* molecules in general. For J PC = 1 ++ , the radiative transitions into 7% c j 
is completely forbidden by C invariance. Even for the molecules of 1 H , the radiative decay 
width should be rather small. Hence a small or negligible radiative decay width poses no 
problem for the interpretation in terms of a molecule. 

If J PC of the most attractive channel of DD coincides with those of a charmonium, a 
substantial mixing can occur between DD and the charmonium [9] . This is a real possibility 

for our most favorite charmonium candidates with J PC = 2 and If X(3872) is such 

a mixed state of a molecule and a charmonium, its radiative decay width is suppressed 
relative to that of a pure charmonium by the fraction of the molecule mixing. 

In fact, this kind of mixing is needed for a very loosely bound molecule state to be 
produced in B decay with a significant rate [10]. The production amplitude of a bound 
state is proportional to the "wave function at origin" \I/ (0) for the same reason that the 
produciton amplitudes of n and K are proportonal to f n and fx, respectively, which are 
interpreted as the wave functons at the origin of quark- ant iquark. For a loosely bound state 
with binding energy AE, |^/(0)| 2 is a small quantity proportional to (mAE) 3 / 2 . For the DD* 
molecule, this is much smaller than |\1/(0)| 2 of the charmonia for which the charm quark mass 
is not sharply defined. If one accepts this argument, production of a pure DD molecule is 
minuscule relative to charmonium produciton in B decay. Physically speaking, D and D 
must fly in parallel with practically zero relative velocities in order to form a molecule in 
B decay. Such a phase space is a tiny, almost negligible fraction of the final-state phase 
space of B decay. The only way to enhance the molecule production in B decay is through 
a substantial mixing with a charmomium. 1 



For production of a pure DD* molecule, a more favorable environment is near the DD* threshold 
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V. DIPION AND ANGULAR DISTRIBUTIONS 



A. Scalar dipion 

The most important information in determining the quantum numbers of X(3872) is the 
spin of the dipion in 7r + 7r~ J/ip. Let us first consider the case that the dipions of / = are 
entirely in s-wave, J p ^ = ++ . If one recalls that the d-wave 7nr cross section is negligibly 
small in this region, this is a reasonable possibility. However, it is experiment that should 
eventually determine whether J^+t,- = or not. If it + it~ is really a scalar dipion, the n + (or 
7r~) momentum in dipion rest frame should show no angular correlation with other vectors; 

dr/<m n ± = -*-r , (8) 

47T 

where fl n ± is the solid angle of p w + — p^ in the dipion rest frame, as measured with respect 
to the direction of any momentum, e.g. the J/ip momentum in the X rest frame. If this test 
shows that the dipion is indeed a scalar, the dipion angular distribution in the rest frame of 
X(3872) happens to be independent of dynamics and unique for J p = l + ~ and 2 since 
the zero-helicity amplitude of J ftp is forbidden for them. We elaborate on this below. 

The X(3872) state is produced in the zero-helicity state in the B rest frame when it is 
produced in B ± — > K ± X(3872). In the X rest frame, X(3872) is in | J, 0) when the quanti- 
zation axis (call it the z-axis) of J is chosen along px in the B rest frame. According to the 
parity constraint of the helicity formalism [11], the helicity amplitudes (s&, A&; s c , A c | J a , M) 
for a — > b + c in the a-rest frame obey 

(s b , X b ; s c , A c | J, 0) = VaVbVc(-^) J+Sb ~ Xb+Sc ~ Xc (s b , -A 6 ; s c , -A c | J, 0), (9) 

where f] a ,b,c are the instrinsic parities of a(= X),b(— 7r + 7r~),c(= J/ip)- In the present case, 
rib = (— lY nn and i] c (= r)j/f) = —1, and s c = Jj/f = 1. We apply this relation to the case 
of a scalar dipion, s& = A& = 0. Denoting the helicity of J/ip by h instead of A c , we obtain 
from the parity constraint 

(0, 0; 1, h\J, 0) = Va (-l) J - h (0, 0; 1, -h\J, 0). (10) 

The factor i] a (— 1) J is equal to —1 for J p = 2~ and 1 + . Therefore the h — amplitudes of 
J/ip must vanish: 

(0,0; 1,0|J,0) = 0, (J PC = 2— and 1 + "). (11) 

Let us suppose that one measures the angular distribution of the dipion momentum = 
Pvr+ + Ptt- in the X rest frame, choosing the z-axis along the X momentum in the B ± rest 
frame. Following the standard formulae [11], we obtain with Eq. (11) 

/ cos 2 0™ sin 2 9 wn , (J PC = 2-) 
rfcos^^lsin 2 ^, (J PC = 1 + -), [LZ) 

where 9 nn is the polar angle of p^. (See Fig.l.) 



in e + e annihilation where the relative motion of D and D is restricted to be small. 
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FIG. 1. The dipion scattering angle 6 wn , which is defined in the X rest frame. 



For the positive parity states other than J PC = 1 + , the longitudinal polarization decay 
amplitude enters the angular distribution: 

dr J 1 (J PC = +_ ), 

dcosfl™ K \ cos 2 9 n7T sin 2 9 n7T + ±k(3 cos 2 9 wn - l) 2 ( J PC = 2 + ~), (13) 

where k = |A | 2 /(|Ai| 2 + |A_i| 2 ) with = \A-i\. The angular distribution for J PC = 2 + ~ 
can mimic that of 2 when k is small. 



B. Tensor dipion 



The zero total helicity state X n7T — h = can be produced as well as the nonzero helicity 
states for the tensor dipion. That is, more than one value of L is allowed for a given J p 
of X(3872) in the case of the tensor dipion. Consequently the angular distributioins and 
correlations are dependent on dynamics. However, kinematics of the X(3872) decay allow 
us to make a special approximation. Since the tensor dipion mass is produced with the 
invariant mass close to the upper end of the kinematical boundary (~ 775 MeV), the dipion 
and J/ip are most often produced nearly at rest in the X rest frame. Let us select those fat 
dipions. We then expect that the relative orbital angular momentum between the dipion 
and J/ip is in s-wave [L = 0) and that the nonrelativistic approximation should be good. 
Therefore, first of all, the angular distribution of the tensor dipion momentum should be 
flat in the X rest frame: 



dT/dcos0™ ~ (l/2)r , (14) 

where the sign of ~ means "in the nonrelativistic approximation" . This may not be an easy 
test since the fat dipions do not move much in the X rest frame so that determination of 
the p nn direction is subject to large errors. 

Eq. (14) only tests L = 0, not directly the quantum numbers of the dipion. We shall 
be able to test the tensor nature of dipion by measuring the angular correlation of p„- = 
Ptt+ — Ptt- °f the dipion rest frame with the X momentum in the B ± rest frame. Since 
angular momenta conservation holds among spins alone for L = 0, the spin components of 
7T+7T- and J/i/j make up \ J,M) = |2,0) of X(3872) as 

|2, 0) x = (1/V2XI2, +1)^11, -1) J/V , - |2, -1)^11, +l>j /v 0, (15) 
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when the spin quantization axis is chosen along the X momentum in the B rest frame. 
The s z = states are missing in the right-hand of Eq. (15), which is characteristic of the 
Clebsch-Gordan composition involving \ j, m) — |1, 0). This leads to the angular distribution 
of Pt,- as 

15 

dT/d cos 9 n ~ — r cos 2 9 n sin 2 9 W , (J PC = 2'-) (16) 

where 9 n is the polar angle of = p„.+ — p^ in the dipion rest frame as measured from the 
X momentum in the B rest frame. Eq. (16) can be obtained by an elementary calculation 
of the nonrelativistic decay amplitude made of three polarizations, e^ kl e l xS lk n e l j^. The test 
will be done most accurately if one selects the dipions of 700 MeV < m nn < 775 MeV. The 
Belle data show many events in this mass region relative to the region below it [1,2]. 

C. Vector dipion 

The spin parity of interest in the case of the vector dipion is J PC = 1 ++ for the molecule. 
The vector dipion most often forms a p meson. Since p and J/ip are nearly at rest in the X 
rest frame, the relative orbital angular momentum between them should be L = just as 
in the case of the tensor dipion. In this circumstance the same simplification occurs. The 
angular distribution of p^ tests L = by Eq. (14). The vector nature of the dipion can be 
tested by the angular correlation, 

dT/dcos9 n ~ ^r sin 2 ^, (J PC = 1++), (17) 

where 9 n is defined in the same way as in Eq. (16). 

VI. SUMMARY 

The most likely candidate for the narrow resonance discovered by the Belle Collabo- 
ration is the l 3 D 2 charmonium state. Depending on the outcome of the radiative decay 
measurement, however, the 2 1 P± charmonium may become a better alternative. Among the 
molecules, the J PC = l +± molecule states of DD + DD* are acceptable as far as quantum 
numbers are concerned and furthermore would easily satisfy the constraints imposed by the 
absence of radiative decay mode. On the other hand, the final states of 7r + 7r~ J/ip with a 
scalar dipion and a vector dipion cannot be produced by a single resonance because of charge 
conjugation invariance no matter how badly isospin symmetry is violated. 

The crucial information leading to determination of the quantum numbers of this reso- 
nance should come from spin of the dipion. The various angular distributions and correla- 
tions that have been presented here will help us in identifying the nature of this resonant 
state conclusively. 
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TABLE I. Unnatural spin-parity and charge-psrity state. L stands for the relative orbital 



angular momentum of 7r + 7r and J/tp. 

J™ I L JPC 

0++ 0,2 1 0+-,l+-,2+- 

2 2~ 

2++ 0, 2 2~ 

1 0+-,l + -,2+-,3 + -,4+- 

1— 1 1++ 

1 0-+,l- + ,2-+,3-+ 
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